When dunes move together, structure of deserts emerges 
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Crescent shaped barchan dunes are highly mobile dunes that are usually presented as a proto- 
typical model of sand dunes. Although they have been theoretically shown to be unstable when 
considered separately, it is well known that they form large assemblies in desert. Collisions of dunes 
have been proposed as a mechanism to redistribute sand between dunes and prevent the formation 
of heavily large dunes, resulting in a stabilizing effect in the context of a dense barchan field. Yet, 
no models are able to explain the spatial structures of dunes observed in deserts. Here, we use an 
agent-based model with elementary rules of sand redistribution during collisions to access the full 
dynamics of very large barchan dune fields. Consequently, stationnary, out of equilibrium states 
emerge. Trigging the dune field density by a sand load/lost ratio, we show that large dune fields 
exhibit two assymtotic regimes: a dilute regime, where sand dune nucleation is needed to maintain 
a dune field, and a dense regime, where dune collisions allow to stabilize the whole dune field. In 
this dense regime, spatial structures form: the dune field is structured in narrow corridors of dunes 
extending in the wind direction, as observed in dense barchan deserts. 



In contrast with the layman's view, not all deserts are 
vast sand seas. Depending on the variability of the lo- 
cal winds, sand dunes can adopt various shapes. When 
viewed from above, they mimic large stars, long linear 
ridges or crescent structures [l], Q ■ The crescent shaped 
dune, called barchan fl], is a prototypical model of 
sand dune dynamics [J] and its properties as an isolated 
object are now well understood |1|,|6|. However, barchans 
are usually found in large dune assembly, counting tens 
of thousands of dunes [3, 0, S]- Barchan fields are ob- 
served in regions where a rocky, non-erodiblc floor is 
blown by a prevalent unidirectional flow and are ubiq- 
uitous on Earth, on Mars and even underseas. The very 
existence of a barchan field is in apparent contradiction 
with the fact that barchan dynamic displays an unstable 
fixed point. Dunes will either grow or shrink if their size 
departs from their equilibrium size, set by the balance 
between sand loss and sand capture. In contrast with 
this unstable behavior, dune size ranges from a few me- 
ters to several hundred of meters within a field. The size 
distribution does not display a lack of small dunes or an 
anomalous number of huge barchans 0, 111 E3l ■ Further- 
more, barchan fields may spatially be structured in nar- 
row corridors, which extend in the wind direction. These 
corridors organize the dune field in stripes of dense (resp. 
diluted) barchan areas, where dunes are smaller (resp. 
larger), while neither the local conditions such as wind 
velocity or granulometry, neither the boundary condi- 
tions differ [7|. Thus, it is commonly thought that dune- 
dune interaction, such as dune collisions, are at play to 
sustain a dune field over longtime, to structure the field 
in corridors and to select the dune size. Field studies 



and underwater experiments have shown that du nes can 
indeed exchange sands during collision events 
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Dune collision led to merge and split mechanisms, which 
can be stabilizing providing that large dunes are regu- 
larly split into smaller ones, as proposed using a mean 
field approach [l3|. Although this idea is a first step, 
it is not enough to fully understand how collisions set 
the structure and affect the stability of a large assembly 
of dunes. Numerical studies implementing phenomeno- 
logical rules of collisions have been runned in order to 
forecast the statistical properties of a large dune field in 
which collisions take place Q . If they succesfully recover 
a size selection, the critical aspect of spatial organization 
of dunes in corridors has never been taken in considera- 
tion yet. It is a general problem in which the large scale 
emergent property comes from the complex, local, inter- 
actions of many objects, what perfectly falls within the 
scope of an agent-based model. 

Here we use an agent-based model, implementing dune 
collisions over the whole dune field, to infer the dynami- 
cal statistical behavior of a large dune field in the limit of 
long observation time, its stationary state and the possi- 
ble emergence of spatial structures. In order to identify 
the physical mechanisms involved in the emergent prop- 
erties of dune fields, this model is restricted to the mini- 
mum ingredients of dunes dynamics and interactions. 
The barchan shape is characterized by a low-slope up- 
wind back and an avalanche face downwind, which is 
framed by two arms pointing in the wind direction (see 
Fig. Q2a-b)). The width, length and height of barchans 
are linearly related to each others such that their morpho- 
logical state can be defined by one parameter only. Sand 
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FIG. 1: Elementary rules for barchan collisions (a) Morphol- 
ogy of a barchan. (b) Parameters used to describe a collision: 
two dunes of size Wi and Wj interact when they are closer 
than do and within a cross section of cry . Each dune is losing 
a volume per unit of time $ and it is fed by sand influx q a . 
(c) The four cases of binary dune collisions. Depending on 
the upward projections s = ^jf , the collision is total (s = 1) 

or lateral (s < 1). Then, the downward projections r — ^ L 
regarding to a merging e m or a splitting e s threshold is taken 
into account. The volume of sand is conserved during a colli- 
sion. Grey levels show the sand redistribution. 



erosion and deposition processes force barchans to move 
downwind. Their propagation velocity, which increases 
with the wind shear stress and decreases with the sed- 
iment influx, is inversely proportional to the dune size. 
The avalanche face acts as a sand trap and barchans can 
propagate over long distances without losing much sand. 
Yet, small sand loss occurs at the tip of the barchan 
arms. Starting with a non-null value, the sand loss in- 
creases very weakly with increasing dune size, and can 
be considered as a constant [l4[. On the other hand, 
the input sand flux is proportional to the dune width. 
As a result, the fixed point (where loss and gain are bal- 
anced) is unstable and an isolated barchan can only grow 
or shrink and eventually, disappear Indeed, below 

a critical size, the barchan loses its avalanche face, turns 
into a dome-like structure and quickly vanishes. 

In our model, dunes are described by their width w 
only and, for the sake of simplicity, are cubic. They prop- 
agate downwind at a speed v: 



a 
w 



(1) 



We assume that dunes lose sand homogeneously along 
their downwind face. We call $ the volume lost per unit 
of time because of wind erosion. Barchans can also grow, 



due to an incoming sand influx per unit of length trans- 
verse to the wind, q a . The volume V of an isolated dune 
will then vary in function of time as: 



dV_ 
~dt 



(2) 



Equation [2] contains the fundamental instability of one 
isolated barchan of unstable equilibrium size w = "I>/<7a- 
If the dune shrinks below the dome size w c , it is removed 
from the field. The model does not conserve the mass. 

One hypothesis to reconcile the unstable behavior of an 
isolated dune (eq. [2|) with the existence of dense barchan 
field is to consider dune collisions. Smaller dunes are 
faster (eq. [T]) and can collide with larger, slower dunes 
what leads to a transfer of mass between dunes [^, [l3| ■ 
Dune collisions led to merge and split mechanism depend- 
ing on the relative size of dunes and their lateral align- 
ment. Numerical studies have shown that those parame- 
ters set the result of a collision [15Hl7| . If the incoming 



dune is very small, it is simply absorbed by the larger, 
slower one. If the dunes are of similar sizes, a redistribu- 
tion of mass occurs, and one (or several) small dunes are 
emitted at the- front while a larger dune is formed at the- 
back. 

In our model, two dunes are in interaction if they are 
closer than a distance do in the wind direction and if 
their width projections overlap as shown in figurc[T|3. Wc 
consider two types of interactions: a distant one through 
emitted sand capture and sand flux screening and a close 
one through collision. The distance do reflects a typical 
distance for the sand flux to get diluted laterally. Let's 
consider two dunes i and j, i being the downwind dune 
and j the upwind one. The size of i is noted Wi and the 
overlapped width er^-. We define the upward projection 
s = — and the downward projection r = — . The lee- 

Wj L J Wi 

ward dune catches a part s of the sand lost by the upwind 
dune. In the same time, the upwind dune screens the lee- 
ward dune on the width Uij from any flux coming upwind 
of i and j. If the upwind dune j is the only one that is 
closer than do to the downwind dune i, the volume of the 
latter varies as dVi/dt = s x $ + q a (u)i — a^) — Note 
that these colian mass exchanges do not affect the aspect 
ratio of the dune (cubic) nor their position. 
Two dunes collide when they overlap following the rules 
shown in Fig. [TJ;. When s = l (perfect overlap), r is the 
size ratio. The two dunes merge if r is smaller than the 
merging threshold £ m , and the new dune gets a volume: 



Vi 



t+At 



(wj) 3 + (w*) 3 . When 



1 and r > e m , 



the 



total sand is redistributed into two new dunes of volume: 



v 4 



t+At _ 



<j? and K t+At 



j) - vij^i <*m* "j — "ij< 

When s < 1 (partial overlap), the sand is redistributed 
into two or three dunes, respectively for r values big- 
ger or smaller than the splitting threshold e a . When 
s < 1 and r ) e Si the two dunes exchange sand the 
same way as when s = 1 and r ^ e m . When s < 1 
and r < e s , the bumping dune i is unaffected while the 
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TABLE I: Parameters of the simulations. 



gets a volume: V k 
Vr At = K') 3 - a ijW . 



bumped dune is split in two dunes. The ejected dune k 
' t+At - anwl while V* +At = {w)f and 
Note that the centers of mass 
of the new dunes are set at the barycentric positions of 
the incoming sand, which may shift the dunes laterally. 
To compensate for sand dune loss, dunes can appear by 
nucleation anywhere in the dune field where there is an 
empty place with a probability per unit time and per 
unit of surface A. Their size is arbitrary set to wq. These 
nucleations are the trace of topographical defects that 
promotes sand deposition 0, Q ■ This choice maximizes 
the effect of noise. Thus any emerging behaviors will be 
robust. 

Looking for stationary behavior of large dune field, the 
dune field boundaries are periodic and the field is long 
compared to the typical distance of dune interactions, 
L > do- Since, dunes move along the wind direction, 
we expect that the perpendicular direction I, does not 
play a major role in dune-dune dynamics. The field is 
initially filled with dunes at random positions homoge- 
neously chosen. Their sizes follow a constant probability 
which is centered on wq and with a minimum cut-off value 
of w c . The numerical method used to compute the as- 
sembly of dunes in this large field is based on synchronous 
algorithm and off-lattice dynamics as in self-propelled 
particles models 
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By their non-trivial kinematics, 
dunes can indeed be considered as self-propelled parti- 
cles, which exchange mass - or momentum - with their 
neighborhood. Values of the different parameters used in 
the simulations are reported in table [U In particular, we 
assume that there is little sand around the dunes, so that 
the ambient influx q a is null. It implies that a lonely dune 
always vanishes. From a simulation to an other, the dune 
density of the field is trigged through a fixed nucleation 
rate A and a changing erosion rate $. 

From the microscopic parameters, seven independent 
dimensionless numbers can be built. One can define three 
length ratios (wo/do, w c /do, w/do = &/(q a do) — infinite 
here), e s and e m which control the dynamics of collisions 
and two times ratio. The three relev ant times are ex- 
pected to be the time of disappearance of an isolated 
dune of size wq: i 00 i = (wq — 1^)/$, the typical nucle- 
ation time: t nuc = (A x Oq) and a collision time: t co \, 
which could be evaluated as the smallest time for two 
interacting dunes to collide i co i = do/ (a/w c — a/wo). In 
the present study, all parameters but $ are kept con- 
stant (see table |T|: t nuc = 2048 and t co \ ~ 11 while 
tool £ [20;6.710 3 ], an isolated dune of size wq travels a 
distance 2 to 660 times its initial size before disappearing. 
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FIG. 2: Characteristics of diluted and dense dune field, (a) 
number of dunes along time for dense and dilute regime, (b) 
size distributions for dilute and dense regime. The continuous 
green line is the analytic law of diluted distribution, eq.2] (c) 
mean width (w) and (d) variance of the width E 2 versus £ for 
different system lengths. Color online. 



We define the control parameter £: 



4> 



\dl 



(3) 



It measures the balance between disappearance of dunes 
due to loss of sand and dune nucleation. So, one can 
expect that a field gets emptied when (Cl. 

Interestingly, the dune field always reaches a station- 
ary state within our range of parameters (see Fig. Wfr)- 
For small £, the barchan field is diluted, a few dunes are 
dispersed across the whole field and dune collisions are 
rare. Dunes can be considered as separate, unstable ob- 
jects whose disappearance is balanced by nucleation only. 
£ is actually the exact dimensionless stationary density 
in the limit of £ <C 1. The normalized distribution of size 
P(w) (see Fig. [2)}) follows the analytic distribution: 



P{w) 



3w' 2 



for w £ [w c ; Wq] 



(4) 



which can be derived from the individual dynamics 
(Eq. [21). Therefore, the typical size of dune is about 
3/4 wo, when w c is small enough. Such a field with low in- 
teraction between dunes compares to diluted deserts such 
as the barchan field of La Pampa de la Joya in Peru 0] ■ 
At large density £ > 1 collisions dispatch sand in a 
non-trivial manner. The field state remains stationary 
(Fig. [2ja,), but the size distribution is fundamentally al- 
tered, and shifted to small sizes (Fig. [2)3). Dune collisions 
tend to increase the number of dunes more rapidly than 
the effect of nucleation itself. With e s and e m set to 
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FIG. 3: Spatial organization of a dune field (a) Density of 
clusters p c \ and probability P n for a dune to be in a cluster vs 
£ for different system sizes, (b) snapshot of a computed desert 
at ^ = 2.9, colors stand for dunes size, (c) closeup of a cluster 
of sub-figure (b). (d) Profiles of mean size (w) and number 
iV of dunes show anticorrelation, see e.g. dashed and dotted- 
dashed lines. Profiles are computed by averaging in the bold 
dashed box shown on figure (b). (Color online. Movies are 
available as supplementary materials) 



0.5, the effective predominant collision type is the frag- 
menting one that creates an additional dune (Fig. [Ip). 
Note that, the most probable size is close but larger than 
the minimal size of dunes w c (Fig. EJd) . A dense assem- 
bly of dunes is not a trivial homogeneous field with fre- 
quent collisions. On the contrary, dense spatial clusters 
of small interacting dunes (i.e. inter-dunes is smaller 
than do) develop and gather the major part of the dunes 
as seen on Fig. [3l These cluster, with sharp boundaries, 
are elongated in the wind direction, so that the field is 
self-structured in a corridor-like pattern where the local 
density is a highly fluctuating quantity. If we restrict our 
measure to local low density, the dune size distribution 
shows a maximum at Wq as in a diluted desert: the lo- 
cal dune size is directly correlated to the local density of 
dunes. A dense field looks like a dilute field of big dunes 
with dense corridors of small dunes (seef Fig. [3ji) . These 
spatial structures and relation between dune size and di- 
lution are similar to what is observed in the long barchan 
field that extends in the Atlantic Sahara (Morocco, see Q 
and supplementary material). Since many small and fast 
dunes are created, another observed effect of collisions is 
the spreading of the speed distribution. They can impact 
larger dunes and lead to a succession of avalanche-like col- 
lisions. Thus, transitory times are very different between 
diluted and dense deserts. Whereas the number of dunes 
relaxes normally in diluted deserts, dense deserts exhibit 



nearly periodic blow-up of their population. 

We identified two stationary states (at £ 1 and 
£ > 1) where the number of dunes but also the size dis- 
tribution, the spatial arrangement in the field and the 
relaxation to the equilibrium are dissimilar (fig. 2 and 
3). This is the signature of two different field dynam- 
ics: a diluted field where dunes barely interact, and a 
dense field whose dynamics is controlled by dune colli- 
sions. It is therefore legitimate to ask whether there is 
any phase transition when transiting from one to another. 
We checked that the mean dune size and its variance 
change continuously when £ was varied (Figl^F and d). 
More generally, whatever the order parameter we looked 
at, we found it to change smoothly without any diverging 
moment. Furthermore, we did not detect any finite size 
effect (see Fig. [5p andOJt), which could sig n a continu- 
ous [19( or a first order phase transition [20j . We neither 
found any influence of initial conditions to the final state 
nor detect any meta-stability. Even if the two stationary 
states are very different, there is no phase transition but 
rather a smooth cross-over when varying £ in this set of 
other fixed parameters. Note that the system we stud- 
ied is far from equilibrium. Thus a phase transition was 
allowed even at low dimension, in contrast with systems 
at equilibrium [2l| . 

In conclusion, wc introduced a minimal agent-based 
model of barchans in interactions, in which kinematics 
and interactions are set in considering experimental ev- 
idence of dune collisions. Its domain of validity cannot 
extend outside the framework of asymptotic limits: in- 
finite size and infinite time of observations. However, 
varying the life time of barchans due to sand loss, we 
showed a smooth cross-over between a diluted desert to 
a dense desert where dunes aggregate in elongated clus- 
ters. Computed deserts self-organize in corridor-like pat- 
terns whith dense regions of small dunes, and diluted 
spaces of larger dunes. This is observed in Earth dense 
barchan fields. Our model, although minimal, was able 
to capture the emergence of such heterogeneous pattern- 
ing. In clusters, the typical barchan size is the result of 
avalanche of collisions. Therefore, we demonstrated that 
a dune fragmentation mechanism, which seems to lack 
in previous studies, is a key process in setting the emer- 
gent properties of barchan dune fields. This fragmenta- 
tion mechanism is here provided by collisions, but one 
could expect that other barchan destabilization mecha- 
nisms could play a similar role. 
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